Do we need Feynman diagrams for higher orders perturbation theory? 
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We compute the two loop correction to the beta lunction for Yang Mills theories in the background 
gauge field method and using the background gauge field as the only source. The calculations are 
based on the separation of the one loop effective potential into zero and positive modes contributions 
and are entirely analytical. No two loop Feynman diagrams are considered in the process. 



PACS numbers: 11.10 Ef, 11.10 Gh, 11.10 Hi, 11.15 Bt. 

■ I. INTRODUCTION 

_ , The instanton approach for SU{N) gauge theories with or without fermions has been initiated by 't Hooft [l| and 

■ further developed in and Q . In this method the separation of quantum degrees of freedom into zero modes (spin 
I i dependent) and positive modes (spin independent) is crucial. Moreover the so called zero modes have an " antiscreening 

■ " effect which is ultimately responsible for asymptotic freedom. The presence of fermions has an opposite effect. In 
^ \ Q we suggest that in essence the magnetic properties of the QCD vacuum play a decisive role in the chiral symmetry 
I— {■ breaking. Furthermore we show in [5| that in the process of gluino decoupling from supersymmetrie QCD separation 

' into zero and positive modes is very important. 

A very useful method for computing beta functions for the gauge coupling constant is the background gauge field 

O I method @ which is based on the decomposition of the gauge field into a background gauge field and a fluctuating 

1^ . field, the quantum gauge field. Even from the dawn of this method the background gauge field was regarded as an 

■ alternate source. However the regular sources J{x) and r]{x), ri'{x) (corresponding to the quantum gauge fields and 
\ ghost respectively) are introduced and one uses the conventional functional formalism to derive beta function or other 
. loop corrections. The reason is simple; the background gauge field does not couple linearly to the other fields (as 

^ ■ linear terms are canceled) and it is not obvious how one can compute simply Green functions with the background 

OO , gauge field as a source. 

^^O ■ In the present work we determine the two loop contribution to the beta function for Yang Mills theories using the 

, background gauge field as the only source present in the functional formalism. Of course the beta function is known 
up the fourth order in the MS scheme so our main interest lies in the method that we introduce and the possibility 

\ for that to be developed for higher orders. We rely on the well-known result of the one loop effective potential (derived 

. either in the perturbative or in the instanton approach) and on the decomposition of the one loop operators into spin 

' dependent and spin independent operators corresponding to each field. Our derivation is entirely based on an analytic 

. . , functional approach (see the Appendix) that does not involve the computation of any two loop Feynman diagram. 



X: 



II. THE METHOD 

The Yang Mills Lagrangian in the background gauge field method (where the gauge field is separated into 5° + A° 
and 5° is the background gauge field) has the expression: 



C = -^[F;^ + D^Al-D,Al + r^^AlAlf- 



- ±^{D^^Alf + -c^[{~D^r ~ D^r^-^AlY (1) 

This lagrangian contains quantum gauge fields and ghosts c", and can be separated into a quadratic contri- 
bution and a higher order one. The procedure for extracting the quadratic terms in this lagrangian is standard and 
after integration leads to the one loop effective potential Q. 
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After some simplification tfie quadratic part of the Lagrangian can be written as: 

1 

The trilinear and quadrilincar contributions are summarized below: 

1 



3,4 = -^iD^Al~D,A-^)r'-A"^A'-' - 

- -^r^T'^'^AlAlA'^^'A'"' + c''{-D^r^''Al)c'' (3) 



Then the effective action reduces to 



= j VAVcexp[i J d'^x{C + Cct] = 
= 12?] exp[* I d'^xi-^iF-^ f + Cct +C2+ £3,4]] = 
= J VAVceM^ j rf'^h^(^;.)' + ^2 + C,t]] 

X [l + i j d^xC^A ^IJ J d''xd^y£3,4{x)C3Aiy) + ••••] (4) 
where Cct is the counterterm Lagrangian. 

The quadratic term alone leads to the one loop effective potential whereas using standard functional procedures one 
can find higher order Feynman diagrams and solve for the higher order contributions to it. Since the first two orders 
of the beta function are renormalization scheme independent one may wonder if it not possible to practically deduce 
two loop contributions using only EqQ and an entirely functional approach without considering and calculating any 
two loop Feynman diagram. In what follows we will show that this is indeed the case by considering the background 
gauge field as a source in the functional approach. For that first we need to express in a suitable form the quadratic 
operators: 

- ^[-{D^r - 2r''^F^'^n = --^[-92 + Ai + A2 + A''] = --^A, (5) 
2.9 2(7^ 2(7^ 

where A"'^ + A^ is a spin independent operator and A'' is a spin dependent operator, 

A^ = i[d''Blr'" + Blr'"'d''] 

= -2/"'"=^''^''. (6) 



(7) 



No spin dependent operator acts on ghosts. 

The one loop effective potential for a Yang Mills theory is obtained by computing, 

exp[*r[i?]] = exp[i J d4x[-^(F;j2][dct(AG,i)]-'/'det(AG,o), 
where Ac,i refers to the gauge fields and Ag,o to the ghost fields. This leads to: 

r[S] = I d'xiFp' + I ln[det(AG,i)] - ln[det(AG.o)]). (8) 

Each operator A has the decomposition from Eq ©with the following calculated contribution to the one loop 
effective potential: 

If 1 f 

2 / d'^x{F'^v)^ — ^ — [^'^-^^'^ T2~] / zero modes contribution for quantum gauge fields 

If 12 AP f 

— — 2 / d'^x{F'^^)'^ — ^ ^tIq^^'^ T2~] / d'^^i^'tlv)'^ positive modes contribution for quantum gauge fields 

4(7 J ^ o J 

If 11 M'^ f 

-—-2 / d^x{FI^^)'^ — ^ --[--iVln — / d'^x{F'^^Y positive modes contribution for ghost fields. (9) 
4(7 J ^ o tz J 
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However it is more convenient for us to represent these results as, 

ln[det(Ai + A2)Ga] = JiVX J d^PP^ 
ln[det(A'^)G,i] = J d''x{Fp^ 

ln[det(Ai + A^)gM = ^NX J d^x{F^^,)\ (10) 

to keep track of the proper regularization procedures. The exact role of X will be revealed in section VI. 

Note that in the case of dimensional regularization each term in Ea P^ will be multiplied by a different X factor 
such that the infinite parts are the same whereas the finite parts are different. 

In the next section we will implement a procedure for computing the two loop beta function by considering the 
background gauge field as a source. Knowing the one loop effective potential before and after the path integration we 
use simple results as those in the Appendix to derive the relevant contributions. The key point in the whole approach 
is the separation of the one loop operators in spin dependent and spin independent ones both in the integrand and in 
the final results. In order to justify that let us consider the general form of the integrals with Dij and Bij the spin 
independent and spin dependent operators respectively: 

/ []de.expK,(i?.,+Aj)e.] = 

I J|dxfcexp[-x,(C'*(B + £i)0),ia;,] = V7r[det[0*(S + i:>)0]]-i/2 ^y^^ 

k 

Here O is the orthogonal operator which realizes the diagonalization. Taking the logarithm of the expression we 
observe that at first order at least there is no interference between the eigenmodes of D and B such that 0*{B + D)0 = 
Bd + Dd + ■■■ where Bd and Dd are the diagonalized operators. This means that at least in the first order the operators 
B and D are diagonalized by the same orthogonal matrix. We will use this feature in our computation. 

III. THE QUADRILINEAR TERM 

We start with the simplest contribution to the effective potential, the quadrilinear term; 

J VA^r'^r'^^A^AtA'^^A^'^ exp[z J d^x[-^iFp^ + C2]] (12) 

The correct structure can be obtained from 

.2 r x2 



J VA\ J d^xd^ySix - y) ^^,J^^,,^^^^ eM J d^x^r^F^^'^A^^Al] 
= J VA^r''T''''AlA''^A'^>'A'"'cxp[J d'^xA^r'"'F^'"'AlAl] 



(13) 



so it is clear that this operator comes only from the spin dependent part in the one loop effective potential. Then 
quite clearly using the results in the appendix one finds: 

J VA^r'^r^-AlAlA'^^A^'^cxvli J d'x[~^iFp^ + C2]] = 

z.g2[_i J d^xd^ySix - y) J)Fb>..^y) exp[-2NX J d''zF\z)]exp[^NX J d^zF\z)] + 
+ j d^:cd^y5(x-y)y^j^^^exp[-2NX J d^zF^{z)] x exp[|^A^X J d^zF'{z)] x 
X- — , expf— 2A^Ar / (i'*zF2(z)] X one loop contribution 

= ig^iX^N'^ J d*zF^{z) X one loop contribution x exp[-^NX J d'^zF^{z)][l + ...] (14) 

For each stage of the calculations we select only the term proportional to the background gauge field invariant 
p^^p'^f^'^ and drop conveniently terms proportional to higher order invariants. 
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IV. THE TRILINEAR PURE GAUGE TERMS 

This contribution corresponds to the term: 

/ d'x{A^''DpAtA"^Tn{x) J d'y{A^^D,AlA^^^f^'^-)iy)] exp[z J d'x[-^^{F;,f + £2]] (15) 

It is simpler in this case to work with the gauge tensor F^^t"" where is the generator in the adjoint representation 
such that, 

Tr{F'i^XF'""'n = NFli^F"'"' (16) 

and Bt" = bt^'t". 

First we notice that the part of the term in Eq ()15p that contains covariant derivatives can be easily obtained from 
the spin independent quadratic operator in accordance to: 

(5[f dSexp[-^A"'(Ai + A2)'^M^]] 1 i r , . . 

''sIb^I^ = ^^^^^P^J^5"'^«-P[-^ J 'i'-Ki^' + ^'rAt] (17) 



Then, 



exph^ / d^xA^li^'+I^.TAl] 



de 



,^AlAl5P-5ad5ce + j^AlDpAlA^D^Al] x 

X exph^ j d^xA^iA^ + A'rAl]. (18) 

We need two more component gauge fields which can be simply obtained from the spin dependent operator. The 
desired result is finally obtained from: 



= [^A"'PA"-A''''Al6P^d''''SadSceS{x -y) + ^{A''''DpAl){x)A'"P{u){A''^D,Al){y)A"^u)] x 
9 9 

X cxp[-^ j d^xAl{A^ + A2)A^ + 1 1 d^xAlF!i:^Al]. (19) 

In order to get the desired term we need to multiply by f'^'^'^fden \(<q^({ to the cancelation by symmetry 
of the unwanted first term in the last line of Eq (|T9|) . 

Eq (|19p has the correct structure except for the space time dependence. We will use a small artifice in order to 
correct that. First we write: 

/ d4a;A™''(a;)F™(x)A"'^(x) j d'^ud\A"'P{u)FJ^^'\u)A'"'{v)S{u ~ v). (20) 

Then, 



g2j pa V . V . g2 



^ d\d'^vA"'P{u)F"^{u)A'"'{v)S{u-v) = \a"^p{wi)A'"'{w2). (21) 



S{Siwi-W2)SF^"{wi)g^J y , pa y , X / V , ^ 

The desired contribution can be computed from: 

d^x{A'^^DpAlA^Pm{x) j dSA'^^^D^AlA-^ f^^){y)]cMi j d^x[~^{F;^f + C^]] 
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-^g^ J d^xd*yd\d\S{x ~ u)S{y - t;)/"^'"/''^" x 

exp[— J d''.TA™''i^™"A"'"] X one loop ghost term. (22) 



(5i^™"(M)(5((5(M - w)) 
Using Ea (P^ and the resuhs from Appendix A we obtain: 



-^[J d^x{A^-'DpA:A"^Tnix) J A(A'*'^i5,^^A"'^/^^")(y)]exp[z J d'x[-^iFp^ + C,]] = 



2 



2 

^ ' d'^xd'^yd'^ud'^vSix - u)S{y - x 



2^ ,5(BP)ac(2:)J(B-)de(y) 

^2 



■exp[yX / d^xF^] X 

, — - X exp[-27VX / d*xF'^] x exp[-iVX / d'^aiF^l x ghost contribution 
SF^^iu)d{d{u~v)) i ' '3 J 

= -Aig^N'^X'^ j d*xF'^{x) x one loop contribution. (23) 

V. TERMS THAT INCLUDE GHOSTS 

There is one quadratic term which contains ghosts and two higher order contributions. We will need to determine 
two terms, respectively: 

-^D^c^r'^Alc%D,Aif'''^fA^PAf-)cM^ J d^x[^^{Fp^ + C2]]. (24) 
We start by analyzing the first term in Ea (P^ . 

Both these expressions contain the ghost fields mixed with quantum gauge fields. For the sake of simplicity we 
write: 

c''i-D^'r'"Aly EE Df.cT'^'Alc'' (25) 

which is true up to a total derivative. Moreover the quadratic term must also be written in a similar manner as: 

c^{~D\,c- = {-D^c-c^) (26) 

We can switch in all these terms the order of the ghost field without problem since we are dealing with the square 
of the trilinear operator. Then the analogy with the previous case is obvious and with exactly the same derivation 
we obtain: 

dHc^DpC^A^-^Prm^) j d^yc'^i?.c^A'-r™)(y)]exp[* j d^ x[~ ^{F^^.f + C2]] = 
= [ d'^xd'^yd'^ud'^vSix - u)6{y - t;)/''^'"/'*™ x 



-g^— -—4-^ -— -expf— A^X / d^xF^] x 

2^^ 5(BP)ac(x)6(B-)de(y) ^42 J ^ 



S 



2 



(5F™"(it)5(^(it- w)) 



— exp[-2ArA: / d'^xF^] x 
W J 



5 1 f 

xexp[(-A^ N)X / d*F^] X one loop spin independent gauge contribution 

3 12 J 

= ig^^N^X^ J d'^xF^ix) X exp[-{^NX J d^F^] x one loop contribution. (27) 
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In this approach the second term in Eq (|24p will give no contribution since it will appear as a product of three 
functional derivatives corresponding to the spin dependent, spin independent and ghost terms in the one loop potential 
and this would lead to a result proportional to a gauge invariant (in the background gauge field) of order higher than 
two. 



VI. CONNECTING THE DOTS 



We add the results from Eq ([T4|) . Eq p3|) and Eq(|27|) to obtain for the second order correction: 

,17 
"6 



ig^-N^X^ / d^x{F;,f (28) 



Here X is just the result of the regularization at one loop. After taking into account all gauge and internal indices 
X amounts to a one loop scalar integral so one can write schematically for the proper loop result: 

d-'xd-'yF^^^Uix - y)F-''''{y) (29) 

The two loop expression then corresponds to: 

d''xdSF::AUU){x~y)F^^^{y) (30) 



where {UU){x — y) is the result of the scalar two loop diagram with two bubbles and two external legs. But this 
regularized is just the square of U regularized at one loop so practically we do not need it. So finally we will take for 
X the expression: 

We multiply by a loop factor i to obtain the second order contribution to the coupling constant: 

From that the known result for the two loop beta function is obtained: 
Here we defined /3(.g) = j^^JI^- 



VII. DISCUSSION 



It is important to know the beta function for non-abelian gauge theories for several reasons. First the one loop 
coefficient of beta function was the main clue that these theories are endowed with asymptotic freedom. Second higher 
order coefficients can reveal information about the phase structure of these type of models. And it is always useful to 
learn more about the mathematical structures that lie at the basis of contemporary particle physics. 

In the present work we do not aim to obtain higher order correction to the beta function but rather to introduce 
a new method that can ease their calculation. All previous computations of order higher than two of beta function 
rely heavily on computer technique, rightly so because for example the fourth order involves about 50.000 Feynman 
diagrams. We determine the two loop coefficient of the beta function as a check of the method. In the process we do 
not compute any two loop diagram and the computation in entirely analytic. 

Our procedure can be easily extended to higher orders and be adjusted to work with dimensional regularization 
and in MS scheme. This would amount to small changes in the general set-up. We leave all these for future work. 
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Appendix A 

First let us review some basic steps regarding integration in the functional approach. We start with the simple 
formula: 

Yl / d^k exp[-$,Sy<fj] = n / ^^f" ^M-bixj] = n V r = const[det B]-^/^ (Al) 
k •' k I * * 

In what follows we will drop the constant factors. 

We extend this to a slightly more complicated case; assume the following: 



n J dik eM^^B,M exp[6A,0] = 11 ^ dct(iJ + D) ^^^^ 

Note that B and D correspond in our case to the spin dependent and spin independent operators respectively in 
the quadratic part of the Lagrangian. Let us now differentiate one of the above factors with respect to a quantity 
Hm where the index includes any type of subscript (Note that in the end this Hm will be the background gauge field 
tensor or any component of it). 

-i[det(S + i?rV2^i|^ 1 (A3) 

2 ^ SHrn d, + hi 

% 

Here di and hi arc the eigenmodes of D and B. The operators B and D are diagonalized together in the one loop 
effective potential such that, 

„2 



-C.AjCj ~d,xi - x^K[jX, (A4) 

where and hi are proportional to the square of the background gauge field tensor [F^^^Y whereas Kij and K^^ 
contain higher powers of the background gauge field. It is perfectly safe to ignore corrections of order Kij or K[y 
We need to compute higher order derivatives of various types: 



r g2 
Jl j d^kexp[-C,B,j^j] ^^ exp[-g,AjCjj 



k 

1 6^d, 1 Sd, 5dj ,1 



2SHnSH„,h,+d, 5Hm5Hn^A{h+di){bj+dj) 2{h,+d,Y 

:[dct(B + D)]-i/2. (A5) 



Furthermore, 



/5 ^2 



S^d, ,1 1 11 



(771— TTTTI— r-TT + oTI— r-7T2<5*fe) + 



'5Hp SHraSHn^A (6, + d,){bj + dj) 2 {h + diY 
Sbk Sd, Sdj 1 1 1 1 1 1 



SHpSH,nSH„' 8 {h, + d,){b, + d,){hk + dk) A {h, + d,)^bj + dj) A {b, + d,)^{bj + d,) 
We treat separately the ghost terms. Thus, 
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from which we can deduce, 



Yl f dOkdel-^ e^p[0*C^J9^] = n / dz^dzl-^ exp[-z*ciz^] = j^^iYlck] 



(AS) 



and further, 



SHmSHn Ci SH,n 5Hn cf CiCj ^ 

It is useful to give in what follows some results regarding the functional derivatives of the square of the gauge tensor 
(Here F^^^f^ = F^^, where is the generator in the adjoint representation). 

' d^zTr{F^,f{z)] = 



5BL{x)5B-,Sy)l 

i[{Fp„)ae{x)5cd5{x - 2/) - {Fp^)dc5ae5{x " v)]- (AlO) 



Furthermore from this one can deduce: 

d'^xd'^yd^d^vSix - u)S{y - v)f^^"^fen ^ 

exp[fcl / d'ziF^^^F^'^ni^)] X 



S{BP)acix)S{B-)M 

(52 



SF"^-^{u)S{S{u- v)) 



— exp[fc2 J d^z{F;:,F^n{-^ 



-mkik2 J d*x{Fp^ X exp[(fci + fc2) J d4^(F;,F"'^'')(z)] (All) 
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